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Regular motions in double bars. 

I. Double-frequency orbits and loops 
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ABSTRACT 

Bars in galaxies are mainly supported by particles trapped around stable periodic orbits. These 
orbits represent oscillatory motion with only one frequency, which is the bar driving fre- 
quency, and miss free oscillations. We show that a similar situation takes place in double 
bars: particles get trapped around parent orbits, which in this case represent oscillatory mo- 
tion with two frequencies of driving by the two bars, and which also lack free oscillations. 
Thus the parent orbits, which constitute the backbone of an oscillating potential of two in- 
dependently rotating bars, are the double-frequency orbits. These orbits do not close in any 
reference frame, but they map onto loops, first introduced by Maciejewski & Sparke (1997). 
Trajectories trapped around the parent double-frequency orbit map onto a set of points con- 
fined within a ring surrounding the loop. 

Key words: methods: analytical — stellar dynamics — galaxies: kinematics and dynamics 
— galaxies: nuclei — galaxies: spiral — galaxies: structure 



1 INTRODUCTION 

Double bars in galaxies, also called nested bars or bars within bars, 
are systems with a nuclear bar embedded within a large-scale outer 
bar. Such systems appear to be common in galaxies: recent surveys 
show that up to 30% of early-type barred galaxies contain them 
(Erwin & Sparke 2002), and they are also present in galaxies of 
types as late as Sbc and Sc (Laine et al. 2002). At least 50 double- 
barred systems are now known (Erwin 2004). The relative orien- 
tation of the two bars is random (Buta & Crocker 1993; Friedli & 
Martinet 1993; Wozniak et al. 1995; Erwin & Sparke 2002; Laine 
et al. 2002), which is expected if the bars were to rotate with differ- 
ent pattern speeds. Inner bars, like large bars, are made of relatively 
old stellar populations, since they remain distinct in near infrared 
(Friedli et al. 1996). 

Orbits in galaxies with two independently rotating bars do not 
conserve the Jacobi integral, and it is a complex dynamical task 
to explain how such systems are sustained. This task has been ap- 
proached with TV-body simulations (e.g. Friedli & Martinet 1993; 
Rautiainen, Salo & Laurikainen 2002; Debattista & Shen 2007; 
Heller, Shlosman & Athanassoula 2007), and through studies of 
trajectories with a selected set of initial conditions (El Zant & 
Shlosman 2002). Nested bars influence each other throughout their 
extent: they can change shapes and accelerate as they rotate through 
each other (Maciejewski & Sparke 2000), and the orbital structure 
of the outer bar can be significantly modified by the inner bar, es- 
pecially near its corotation. In general, the piling up of resonances 



E-mail : witold @ astro, ox . ac . uk 



created by each bar may lead to considerable chaotic zones, and, if 
the fraction of chaotic orbits in double bars were large, such sys- 
tems should not last for long time periods. Long-lasting double bars 
should therefore have a substantial fraction of regular orbits, and 
one should be able to describe their dynamics by means similar to 
those invoked in the case of single bars, i.e. in terms of the regular 
orbits that they admit. Such an approach was originally proposed 
by Maciejewski & Sparke (1997, 2000). In this paper we propose 
that the backbone of the oscillating potential of double bars is made 
of double-frequency orbits, and we explicitly show that these orbits 
map onto the loops invented by Maciejewski & Sparke. 

In Section 2, we outline in more detail the approach proposed 
in this paper. In Section 3, we use the epicyclic approximation to 
show that two frequencies on an orbit in the plane of an oscillating 
potential of a doubly barred galaxy are indispensable. We also show 
there that double-frequency orbits map onto loops in the epicyclic 
approximation. In Section 4, we analyze particle trajectories in dou- 
ble bars in the general non-linear case. We develop a method of 
analysis in Fourier space that allows us to distinguish regular tra- 
jectories from chaotic ones, and to recover the three fundamental 
frequencies of the regular orbits. We find that in the general, non- 
linear case double-frequency orbits map onto loops as well, and 
that they can trap around themselves regular orbits. 



2 CONCEPTUAL OUTLINE 

Particle motion in a potential of double bars belongs to the gen- 
eral problem of motion in an oscillating potential (Louis & Ger- 
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hard 1989; Sridhar 1989). In this particular case, the potential os- 
cillates with the relative period of the bars. The best known ex- 
ample of the general problem of motion in an oscillating potential 
is the restricted elliptical 3-body problem. Families of closed pe- 
riodic orbits have been found in this last problem, where the test 
particle moves in the potential of a binary star with components on 
elliptical orbits (e.g. Broucke 1969). Such families, however, are 
parametrized by values that also characterize the potential (i.e. el- 
lipticity of the stellar orbit and the mass ratio of the stars), and the 
orbital periods there are commensurate with the oscillation period 
of the potential. For a given potential, these families are reduced to 
single orbits separated in phase-space. The solution for double bars 
is formally identical, and there an orbit can close only when the 
orbital period is commensurate with the relative period of the bars. 
Howard (1994) found such periodic orbits in double bars with pe- 
riods equal to a multiplicity of the bars' relative period. Such orbits 
are separated in phase-space and therefore families of closed peri- 
odic orbits are not sufficient to provide orbital support for nested 
bars. 

Binney & Spergel (1982, 1984) developed a formalism that 
classifies the orbits based on the frequencies at which they oscil- 
late. The Fourier transform of any dynamical variable on a regular 
orbit (position, velocity) should return a line spectrum. Regular or- 
bits in the plane of a stationary galaxy correspond to oscillations 
with two frequencies: the angular frequency, and the epicyclic fre- 
quency - the orbit's fundamental frequency. Closed periodic orbits 
are a particular case of regular orbits. They oscillate with one fre- 
quency only, whose inverse is the period of the orbit. In general, 
a single frequency in the system may appear when the two fre- 
quencies above are commensurate, or when oscillations with one 
of the frequencies are missing. A rotating asymmetry {e.g. a bar) in 
the system constitutes a periodic driving force. As in the case of a 
driven harmonic oscillator, the orbits always respond to this driving 
force by oscillations with the frequency of the driver. These oscil- 
lations are always present, and if the orbit is to close in the frame 
rotating with the bar, they must be the only oscillations, i.e. the free 
oscillations with the epicyclic frequency must be absent. 

The oscillating potential of double bars changes periodically 
with the relative period of the bars. Since a Hamiltonian system in 
N dimensions with a periodically varying potential is equivalent 
to an autonomous system in N + 1 dimensions (see e.g. Lichten- 
berg & Lieberman 1992, Louis & Gerhard 1988), regular motions 
in the plane of a doubly barred galaxy should correspond to os- 
cillations with three frequencies. As before, Fourier analysis of a 
closed periodic orbit will return only one frequency. Aside from 
oscillations with three frequencies, and with one frequency only, 
there remain oscillations with two frequencies. Since two bars in 
a doubly barred galaxy constitute two periodic driving forces, os- 
cillations with two frequencies should characterize response of the 
system to this periodic driving, when free oscillations are absent. 
Thus, the role of double-frequency orbits in double bars should be 
equivalent to the role of closed periodic orbits in a single bar. In 
the unchanging potential of a single bar, regular orbits are trapped 
around stable closed periodic orbits, which constitute the backbone 
of that potential. Here we postulate that stable double-frequency 
orbits constitute the backbone of an oscillating potential of double 
bars, and trap around themselves regular orbits with the maximum 
number of frequencies, which, in this case, is three. 

The approach outlined above indicates that oscillating poten- 
tials can be supported by orbits that do not have to be closed and 
periodic. With this approach we overcome the problem of scarcity 
of closed periodic orbits that hampered previous searches for sup- 



port of oscillating potentials. The concept of closed periodic orbit is 
too limiting for the investigation of oscillating systems, and we re- 
place it by another, broader description of particle motion, which 
still captures its regularities. Our description, based on double- 
frequency orbits which are not closed and periodic, can provide the 
backbone of an oscillating system for the same reason for which 
closed periodic orbits in stationary potentials do, and account for 
longevity of double bars. 

Double-frequency orbits are difficult to visualize, because 
they do not close in any reference frame. However, in this paper 
we will show that they can be mapped onto closed curves, which 
are much easier to visualize. The mapping is done by recording 
the position of a particle on a double-frequency orbit at every con- 
secutive alignment of the bars. Thus particle's position is written 
only when the oscillating galaxy periodically regains one specified 
shape. In this paper we show that positions of particles on a double- 
frequency orbit recorded at every consecutive alignment of the bars 
populate a closed curve. By construction, this curve transforms into 
itself at each consecutive alignment of the bars. Thus a set of par- 
ticles populating such a curve oscillates periodically with the os- 
cillating potential. Maciejewski & Sparke (1997) discovered such 
curves in a potential with a double bar, and they called them loops. 
Loops can be used as a convenient descriptor of orbital structure 
in an oscillating potential. Like closed periodic orbits in stationary 
potentials, they can indicate the zones that can be supported by reg- 
ular motions. In the remaining sections we will show that double- 
frequency orbits can form the backbone of oscillating potentials, 
and that loops make a unique tool for studying such orbits. 



3 THE EPICYCLIC SOLUTION 

3.1 The general epicyclic solution for any number of bars 

If a galaxy has a bar that rotates with a constant pattern speed, it 
is convenient to study particle orbits in the reference frame rotat- 
ing with the bar. If two or more bars are present and each rotates 
with its own pattern speed, there is no reference frame in which 
the potential remains unchanged. In order to point out formal simi- 
larities in solutions for one and many bars, we solve the linearized 
equations of motion in the inertial frame. As we show below, this 
is equivalent to the solution in any rotating frame and the transfor- 
mation is particularly simple: in the rotating frame the centrifugal 
and Coriolis terms are equivalent to the Doppler shift of the angular 
velocity. 

In a frame rotating with a constant pattern speed Qb , the equa- 
tion of motion of a particle in a potential <3> can be written as 

r' = -V$ - 2(ft B x r') - fl B x (fi B x r'), 

where the position of the particle in the rotating frame is indexed 
by prime. In polar coordinates in the plane of the galactic disc, its 
R' and ip' components are 

R' - JJV + Q B ) 2 = -||, (1) 
R'p> + 2R>(v' + n B ) = -^f^- (2) 

Since the angular velocity ip' in the rotating frame is related to the 
angular velocity tp in the inertial frame by 

ip = ip' + Q B , 

and radial coordinates in both frames are equal (R = R'), equa- 
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R - Rip = 
R(p + 2Rip = 



(3) 
(4) 



tions (1) and (2) in the inertial frame take the form 

.2 _ _<9£ 
OR' 

~Rdcp' 

which are the R and ip components of the equation of motion in the 
inertial frame: 

t = -V<E>. 

This equivalence remains true when vertical motion is included in 
three dimensions. For the rest of this section we adopt the inertial 
frame, in which the R and ip components of the equation of motion 
take a particularly simple form of (|3][4]l. 

We want to linearize equations ([3} and 101 for departures from 
circular motion in a nearly axisymmetric potential. We thus per- 
form two expansions simultaneously. One is a linearization of the 
quantities describing the motion by their expansion to first-order 
terms in R — Ro, where Rq is the radius of the circular orbit in the 
axisymmetric potential. The other is expansion of the potential in 
a Fourier series, with the zeroth-order term being the axisymmet- 
ric component. In both expansions, zeroth-order terms are indexed 
with and first-order terms are indexed with /. These expansions 
can be written as 



R(t) = R + Ri(t), 
ip(t) = ipoo + Hot + ipi(i), 
$(R,tp,t) = <f> (R)+&i(R,<p,t), 



(5) 
(6) 
(7) 



where second- and higher-order terms are neglected. In the iner- 
tial frame, the angle coordinate ip of a particle in the zeroth ap- 
proximation is ipo = ipoo + Hot, where ipoo is the value of ipo at 
t = 0, and where the angular velocity fio on the circular orbit of 
radius Rq relates to the potential $0 through the zeroth order of {3): 
Hq = {l/R )(d$ /dR) \r . The zeroth order of © is identically 
equal to zero. In the zeroth order approximation, the axisymmetric 
part of the potential, $0, is time-independent. The first order term 
<&i, describing asymmetry in the potential, may depend on time. 

In the expansion defined above, the first order parts of {3j and 
© take respectively the forms 



Rj - 4AH Ri - 2R H <pi 
Ro<pi + 2H Ri 



' Q R |fio.?o. 

1 d$i , 

Ro dip 



Ro,vo 1 



(8) 
(9) 



where A is the Oort constant defined by AAHo = J7q — g^r\R ■ 
Furthermore, we assume that, to the first order, the departure 
of the barred potential from axial symmetry can be described by 
a term cos(2ip). If N bars, indexed by i, rotate independently as 
solid bodies with angular velocities Hi, the time-dependent first- 
order correction $j to the potential can be written as 



$i(R, tp, t) = ^ cos[2(v? - Slit)], 



(10) 



where the radial dependence ^i(R) has been separated from the 
angular dependence. The lack of phase in the above trigonometric 
functions means that we define as t — the time when all the bars 
are aligned. Derivatives of J 1 Ob enter the right-hand sides of {8} and 
l[9}, which after introducing 



LUi = 2(f2o — Hi), 
take the form 



(11) 



N 

-R I +4AH R I +2RoH tp I = ^ \ R<) cos{uit+2tp 00 )(12) 



R (p'i + 2Q.OR,! = *i(-Ro) sin(tj l t + 2tp 



(13) 



In order to solve the set of equations d 1 21 1 3b . one can inte- 
grate J 1 3 b and get an expression for Ra<pi, which furthermore can 
be substituted into J12b . This substitution eliminates ipi, and one 
gets a single second order equation for Ri, which can be written 
schematically as 



Ri + KqRi = Mi cos(wit + 29300) + C, f 



(14) 



where Mi 



89 

an \r 



,Kq = 4fi (f2o -A), and C V /2H 
is the integration constant that appears after integrating d 1 3 b - Below, 
we will show that C v can be removed by redefining the guiding ra- 
dius. What is left of J 1 41 . is the equation of a harmonic oscillator 
with multiple driving terms. Thus in the linear approximation, the 
orbital motion of a particle in the potential of rotating bars is equiv- 
alent to the motion of a driven harmonic oscillator. The general 
solution of d!4t is well known and can be written as 

JV 

Ri{t) = Cicos(K t+S)+^2P i cos(ujit + 2<p 00 ) + C ip /K,l,(15) 

i 

where the coefficients Pi are functions of Mi. The first term in the 
right hand side of J 1 5b is the solution of the homogeneous part of 
d 1 4b . and it corresponds to a free oscillation with the local epicyclic 
frequency kq around the circular orbit in the axisymmetric poten- 
tial. The coefficient Ci is unconstrained. The terms in the sum over 
i in J 1 5b are the particular solutions of the inhomogeneous equa- 
tion d 14b . and they describe oscillations driven by the bars. Note 
that the frequency of the driving force (shifted to the frame of the 
rotating guiding centre) is always present in the oscillations. Thus, 
orbits in potentials with N independently rotating bars will oscil- 
late with at least N frequencies, unless the shifted frequencies are 
commensurate. 

The formula for (pi (t) can be obtained by substituting d 1 5 b 
into the time-integrated d 1 3 b . As a result, one gets 



<fii = C2 cos(ft i + 5) + Qi cos(ujit + 2<p o) 



K^Ro 



,(16) 



where again C2 is unconstrained and the Qi are determined by the 
coefficients of the equations above. Note that the integration con- 
stant C^/kq in J15b can be incorporated into Ro in lO, giving the 
new guiding radius Ro + C v /k^. Since A = — |-Ro^§f \Ro, ex- 
panding SI to the first order gives the angular velocity at this new 
radius 



Ho 



Rq + C v /k 



= Ho\r„ H j 



dHo 
~~dR 



\R 



= Ho\ 



2AC V 
k\Ro 



2 AC 

Thus the integration constant j— — in d 1 6b represents the correc- 

K g R 

tion to the zeroth order angular velocity Ho incurred by a change 
of the guiding radius Rq by C v /k,q, which is the integration con- 
stant in d 1 5 b - Hence both integration constants can be incorporated 
by modifying the guiding radius in the epicyclic approximation. 
Eventually, in the epicyclic approximation, the most general solu- 
tion for particle motion in the symmetry plane of a potential with 
iV independently rotating bars is 
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Ri(t) = Ci cos(«;ot + 5) + ^ p i cos(ujit + 2ip 00 ), (17) 



C2 Q 

<pi(t) = — sin(ftoi + <5) + } — - s'm(ujit + 2(poo) +const, (18) 

i 

where Ri and ipi are polar coordinates in the inertial frame. The 
integration constant in d 1 8b is an unconstrained parameter. From 
d!7b and d!8b one can see that up to N + 1 independent frequencies 
are present in this motion. 



3.2 Closed periodic orbits in a single bar 

In a potential with a single bar there is only one term in the sums 
over i in d!7b and d 1 8b . Therefore in general two frequencies are 
present in orbital motion, as expected for motion in a symmetry 
plane of a static potential. However, closed periodic orbits have 
only one frequency. Oscillations with the frequency of the rotating 
bar (corrected for the angular frequency of the guiding centre) are 
always present in orbital motion, analogously to oscillations with 
the frequency of the driving force in a driven harmonic oscillator. 
The linear approximation breaks down when the frequency of free 
oscillations becomes commensurate with the frequency of the driv- 
ing force. Therefore, in this framework, in order for the orbits to 
close, the amplitude of free oscillations with frequency ko has to 
be zero. For Ci = C2 = and N = 1, equations d!7b and l |18l > de- 
scribe closed periodic orbits in the linearized problem of a particle 
motion in a single bar. 

Consider, for the case where free oscillations are absent, the 
change in the values of Ri and ipi after the particle's guiding centre 
returns to its starting point in the frame corotating with the bar. 
This happens after time period 2n/(Qo — fii). Replacing t by t + 
2n/(yi<> — fii) leads to 



fi 



mill — n^2 



(19) 



Ri — P\ cos[u)i(t + ■ 



2tt 



fio — fii 
Pi cos(ojit + A-k + 2<poo 



+ 2<Ax>] 



Thus, after time 27r/(fio — fii) the solution for Ri returns to its 
starting value. The same holds true for tpi and the orbit closes in the 
bar frame. Thus, 2n/(Q,o ~ fii ) is the full period of this particle. 



3.3 Closed periodic orbits in double bars 

While the potential of a single bar is stationary in the bar's frame, 
that of a galaxy with two independently rotating bars is not station- 
ary in any reference frame. For two bars there are two terms in the 
sums over i in d!7b and J 1 3b . and therefore even if the amplitude 
of free oscillations with the frequency ko is zero, two frequencies 
are present in the orbital motion. They are the frequencies of the 
driving forces imposed by the two tumbling bars, and oscillations 
with these two frequencies have always non-zero amplitude. Orbits 
with two frequencies do not close, therefore closed periodic orbits 
should not be the fundamental orbits in double bars. 

However, closed periodic orbits do exist in double bars, and 
in the linear approximation this happens when the frequencies in 
the sums over i in d!7b and d!8b are commensurate, i.e. when 
nuj\ = mu>2 for n, m integer. In that case there is effectively only 
one independent frequency and the orbit closes. However, this hap- 
pens only for a discrete set of values of the angular velocity of the 
guiding centre 



Therefore in the epicyclic approximation there is an infinity of 
closed periodic orbits in double bars, but these orbits do not form 
families continuous with the radius of the guiding centre. Although 
formally any orbit in double bars, which responds to two driving 
forces by oscillating with two frequencies, can be approximated 
by a closed periodic orbit of angular frequency l |19l >, such an orbit 
closes after very many bar periods, and describing motions in dou- 
ble bars in terms of such closed periodic orbits is not very conve- 
nient. Moreover, although any orbit oscillating with two frequen- 
cies has a closed periodic orbit in its infinitesimal vicinity with 
sufficiently large m and n, such a closed periodic orbit cannot be 
considered to be a parent orbit, because orbits with two indepen- 
dent frequencies are not constituted by free oscillations around that 
orbit. Finally, in the non-linear regime, Howard (1994) found that 
closed periodic orbits in doubly barred potentials, with periods pre- 
dicted by < | 1 9b are often unstable, and therefore they cannot serve 
as parent orbits. Thus, closed periodic orbits are unlikely to pro- 
vide sufficient backbone supporting double bars. Below we present 
a more convenient way serving this purpose. 



3.4 Double-frequency orbits in double bars. Loops 

In Sect. 3. 2 we showed that in the case of a single bar, when there 
are no free oscillations, the linear solution d!7b and d 1 8b returns 
only one oscillation frequency, which indicates a closed periodic 
orbit. In a solution for double bars, two independent frequencies are 
present in general, hence one should expect that double-frequency 
orbits take place of closed periodic orbits there. Such orbits do not 
close in any reference frame, because when a term from one bar in 
d!7b and d!8b returns to its initial value, the term from the other bar 
does not. 

However, consider the change in value of Ri and ipi after time 
tv/(Q,2 — fii), which is the relative period of the bars. Every such 
time interval, the mass distribution in the galaxy and its potential 
are the same in the frame rotating with either bar. With amplitudes 
of free oscillations in d!7b set to zero one gets 



Ri = Picos[wi(t + 

+ P2COs[i02(t + 



Pi cos(o;it + 2tt 
P2 cos(oj2i + 2n 



o 2 - Hi ' 
n ) 

Q.2 - fii 



Q.2 

fio 



fii 

0.2 



f 2<poo] 
f 2v?oo] 
+ 2^oo) 
+ 2ipoo) 



2 - fii 

= Pi cos(o;it + 2tt + 2(poi) + P2 cos(tJ2t + 2<poi), 

where (poi = (£00 + 7r ^~^ . The same result can be obtained for 
ifii. This means that the time transformation t — * t + n/(0 2 — fii) 
is equivalent to the change in the starting position angle of a parti- 
cle from ipoo to ipoi ■ Consider the motion of a set of particles that 
have the same guiding radius Ro, but start at various position an- 
gles tpoo . This is a one-parameter set, therefore in the disc plane it is 
represented by a curve, and because of continuity of J 1 7b and d!8b 
and their periodicity at a fixed t, this curve is closed. After time 
n/(02 — fii), a particle starting at angle ipoo will take the place of 
the particle which started at 1^01 , a particle starting at <poi will take 
the place of another particle from this curve and so on. The whole 
curve will regain its shape and position every 7r/(fi2 — fii) time 
interval, although positions of particles on the curve will shift. This 



Regular motions in double bars. I. Double-frequency orbits and loops 5 




Figure 1. Two example trajectories (one in the two left panels, one in the 
two right ones) of a particle that moves in the potential of two indepen- 
dently rotating bars from Model 1 of Maciejewski & Sparke (2000). Large 
dots mark the starting position of the particle at 1 .7-kpc distance from the 
centre of the galaxy, on the minor axis of the bars, at the moment when the 
corresponding axes of both bars overlap. The initial particle's velocity vec- 
tor is in both examples perpendicular to the minor axis, and its value equal 
to v x = 220 km s — 1 (left panels) and v x = 150 km s -1 (right panels). 
The particle is followed for 10 relative periods of the bars, and its trajectory 
is displayed in the frame corotating with the big, outer bar (top panels); and 
with the small, inner bar (bottom panels). Each bar is outlined in its own 
reference frame by a dotted line. 

curve is the epicyclic approximation to the loop introduced by Ma- 
ciejewski & Sparke (1997, 2000): a curve made of particles moving 
in a given potential, such that the curve returns periodically to its 
original shape and position. In the case of two bars, the period is 
the relative period of the bars, and in the epicyclic approximation 
the loop is made out of particles having the same guiding radius 
Ro. Particles on the loop respond to the periodic driving force from 
the two bars, but they lack any free oscillations. Note that the ap- 
pearance of the loops is the same, no matter which bar's reference 
frame is used. Thus, loops are a convenient descriptor of orbital 
structure in an oscillating potential. An example of a set of loops in 
a doubly barred galaxy in the epicyclic approximation was given by 
Maciejewski & Sparke (1997). The loops occupy there a significant 
part of the disc. 

Note that the linear solution presented here gives a prescrip- 
tion for constructing a loop in a general case as well. If a particle 
is on a double-frequency orbit, then the loop onto which the orbit 
maps can be constructed by registering the position of the particle 
every relative period of the bars, until the whole 2n range of an- 
gular coordinates is densely populated. However, this construction 
requires following the particle for many relative periods of the bars, 
hence it will likely recover loops that are maps of stable double- 
frequency orbits only. 



4 FULL NONLINEAR SOLUTION FOR PARTICLE 
MOTION IN NESTED BARS 

4.1 Three types of trajectories 

Tools and concepts useful in the search for ordered motions in dou- 
ble bars are best introduced through the inspection of particle tra- 
jectories in such systems. We use the potential of Model 1 defined 
in Maciejewski & Sparke (2000), so that our results can be directly 
compared to the results there. In order to describe gravitational po- 
tential in their models, Maciejewski & Sparke use the analytical 
formulae originally proposed by Athanassoula (1992). The formu- 
lae comprise the modified Hubble profile representing the bulge, 
the Kuzmin-Toomre profile representing the disc, and the Ferrers 




Figure 2. Most regular trajectory that starts at 1. 7-kpc distance from the 
centre of the doubly barred potential from Model 1 of Maciejewski & 
Sparke (2000), followed for 20 relative periods of the bars, and written in 
the frame corotating with the big bar (left), and with the small bar (right). 
The major axis of each bar in its own reference frame is marked by the 
dashed line. Dots mark positions of the particle at every alignment of the 
bars (see Sect. 4. 4). 

formula describing each bar. The exponent in the Ferrers formula 
is 2. In Model 1 of Maciejewski & Sparke, the outer, main bar is 
7 kpc long, and the inner bar is 60% in size of the outer one. The 
axial ratio of each bar is 2.5. The quadrupole moment of the outer 
bar is only 2.25 x 10 10 M Q kpc 2 in Model 1, and the mass of the 
inner bar is 60% of that of the outer one, hence both bars in this 
model are rather weak. Pattern speeds of the bars are not commen- 
surate (fii = 24.03 km s _1 kpc -1 , fi 2 = 41.90 km s _1 kpc -1 ), 
and indicate rapidly rotating bars: the corotation radius of each bar 
is about 20% larger than its semi-major axis. The Inner Lindblad 
Resonance of the outer bar is located at 2.9 kpc, hence there is no 
resonant coupling between the bars in this model. 

Consider a particle moving in this potential inside the corota- 
tion of the inner bar. Simple experiments exploring a limited range 
of initial conditions - with the particle starting on the minor axis 
of the bars at the moment when the axes overlap and with velocity 
vector perpendicular to that axis - show that if the initial velocity 
is small enough, the particle usually remains bound. Its trajectory 
may look regular or irregular. To illustrate this, in Figures 1 and 
2 we present three trajectories of a particle starting at 1. 7-kpc dis- 
tance from the centre of the galaxy with three different initial ve- 
locities. A typical irregular trajectory is shown in the left panels 
of Fig.l (starting velocity v x = 220 km s -1 ). Since the trajectory 
depends on the reference frame, it is plotted twice, once for the 
reference frame of each bar. However, particle trajectories usually 
look more regular: an example is given in the right panels of Fig. 1 
(starting velocity v x = 150 km s" 1 ). These trajectories look like if 
they were trapped around some stable regular orbit. 

Fine adjustments of the initial velocity lead to a highly har- 
monious trajectory (Fig.2, starting velocity v x = 180 km s _1 ), 
which appears even more regular than the trajectory from the right 
panels of Fig.l. Moreover, this is the most regular trajectory that 
we can find by varying the initial velocity. In particular, we do not 
find any closed periodic orbits that this trajectory can be related to. 
Below we will formalize the meaning of being 'most regular' in 
this context. The trajectory from Fig.2 looks like a loop orbit in a 
potential of a single baiQ (see e.g. fig. 3. 7a in Binney & Tremaine 
1987). In Sect.4.2, we will show that this morphological similarity 

1 in the existing nomenclature, which we continue to use, a loop orbit is 
an orbit in a time-independent potential, and it has no relation to a loop, 
which is a map of a double-frequency orbit in a time-dependent oscillating 
potential 
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angular frequency 

Figure 3. One-sided power spectral density (power spectra) from fast 
Fourier transform of all six dynamical variables described in the text on 
an exemplary regular trajectory. The vertical axis is logarithmic, with the 
value increasing by one order of magnitude every large tick mark. For clar- 
ity, the spectra are shifted in power from their neighbours by 7 large tick 
marks each. 

occurs because in both cases orbital oscillations with two frequen- 
cies are involved. The loop orbit in a single bar forms when a parti- 
cle oscillates around a closed periodic orbit. Therefore, it involves 
two frequencies: the frequency of the free oscillation, and the fre- 
quency of driving by the bar. In Sect.4.2. 1 we will demonstrate that 
in the trajectory from Fig.2, the two frequencies are the frequencies 
of driving by the two bars, while the free oscillation is absent. This 
is the reason why the trajectory from Fig.2 appears most regular: 
adding the free oscillation with its own frequency will make the 
appearance of the trajectory less regular. 



4.2 Orbital frequencies 

In order to see what frequencies are present in the motion of a 
particle on a given trajectory, we Fourier-analyse any dynamical 
variable on this trajectory. Fast Fourier transform method has been 
employed, and the trajectory was traced for 100 alignments of the 
bars, with dynamical variables sampled N = 2 13 times at equal 
time intervals. We transformed the polar coordinates of the par- 
ticle, as well as its polar velocity components. The values of the 
radial velocity v r entered the transformation directly, while the ra- 
dial coordinate r had its average over N samplings subtracted be- 
fore it was transformed. Similarly, the average over N samplings 
of the tangential velocity vt was subtracted from its values before 
they entered the Fourier transform. In order to transform the tan- 



gential coordinate tp, we also had to subtract the component re- 
flecting its monotonic growth. In effect the following values were 
transformed: v r , r - (Er*) /N, v t - (Euj)/iV, tp c - (E<pi)/N, 
where ip c = tp — Qt, t is the time since the moment when the par- 
ticle started, Cl — {tp N — tp 1 )/T, and T is the total time for which 
the particle has been traced. The sum E is over all N samplings, 
and the upper index counts the samplings. Moreover, the angular 
coordinate tp and the tangential velocity v t were transformed for 
the frame of reference of each bar (tpd, tp c2 and vti, Vt%), hence 6 
dynamical variables were transformed in total. 

In Fig.3, we present one-sided power spectral density (for 
brevity called here power spectra) from fast Fourier transform for 
each variable on an exemplary regular trajectory. They are essen- 
tially line spectra, as most of the power resides in a few separated 
peaks. Note the quality of the spectra: the power density in the 
peaks is several orders of magnitude higher than that in the base- 
line. Power spectra of all six dynamical variables transformed have 
the same distribution and similar strengths of major peaks. Only 
spectra of the angular coordinate tp are of slightly lower quality, 
with higher baseline, and power increasing towards the lowest fre- 
quencies. This is likely the effect of inaccuracies in double cor- 
rections applied to the measured value of those variables. Even in 
these spectra, however, the peaks are at least two orders of mag- 
nitude in power over the highest values of the baseline. In the rest 
of this paper, we focus on the analysis of the spectrum of the ra- 
dial coordinate only, for which our Fourier algorithm returns data 
of best quality. 

4.2.1 Most regular trajectory as a double-frequency orbit 

In Fig. 4, we present power spectra for the three trajectories from 
Figures 1 and 2. We display the power density in both the linear 
scale (left panels), and the logarithmic scale (right panels). The 
power spectrum for the most regular trajectory from Fig.2 is shown 
in the top row. It shows sharp peaks at two frequencies: UJ2 = 163.7 
kms -1 kpc -1 and^i = 199.4kms -1 kpc -1 . These peaks contain 
about 95% of total spectral power. Most of the remaining power is 
in secondary peaks at frequencies which are functions of wi and 
u)2, namely wi + u>2, 2a>i, 2^2, W2 — ui etc. The power in the sec- 
ondary peaks rapidly decreases as a function of their frequency: the 
group of three peaks between 300 and 400 km s -1 kpc -1 contains 
less than 2% of the total power, and further peaks contain still one 
order of magnitude less of the total power. 

The two frequencies, u)\ and u>2, at which almost all of the 
power resides, are separated by 35.75 km s -1 kpc . This is ex- 
actly twice fl2 — fli, the difference between the pattern speeds of 
the bars in Model 1. In Sect. 3, we showed that in the linear ap- 
proximation the solution for the radial coordinate dl7t of a particle 
moving in the potential of two independently rotating bars, in the 
absence of free oscillations, consists of two oscillations with fre- 
quencies uj\ and u)2 defined by i ll It . The difference between those 
frequencies is also 2(fi2 — fii). Thus, we imply that the frequen- 
cies of the highest peaks in the power spectrum of the trajectory 
from Fig.2 are in fact the frequencies of driving by the two bars, 
and that they are nonlinear equivalents of the Ui's defined by i ll It . 
From i ll It one can obtain the underlying orbital frequency of the 
particle: Q — 123.75 km s _1 kpc -1 . It is close to its epicyclic ap- 
proximation at the guiding radius ro = (Er l )/jV = 1.803 kpc, 
which is 121.5 km s -1 kpc -1 . This consistency further confirms 
our finding. Thus in both the linear case (epicyclic approximation) 
and in the general case, we constructed double-frequency orbits in 
an oscillating potential of a double bar, with frequencies equal to 
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Figure 4. Power spectra of the radial coordinate on three representative trajectories, shown in linear (left-hand panels) and logarithmic (right-hand panels) 
scale. The plots in the top row show the power spectrum of the double-frequency orbit from Fig. 2. The plots in the middle and bottom rows show, respectively, 
power spectra of the trajectories from the right-hand and left-hand panels of Fig. 1. The particle was followed for 100 relative periods of the bars, with dynamical 
variables sampled N = 2 13 times at equal time intervals. In the bottom-left panel, the thick line marks an additional power spectrum, for the same initial 
conditions, but with the particle followed for 600 relative periods of the bars, with dynamical variables sampled TV = 2 16 times. 



the frequencies of driving by the bars. Note that since the Fourier 
analysis requires the particle to be followed for many relative peri- 
ods of the bars, and since no fine adjustments of initial conditions 
were made, this double-frequency orbit is most likely stable. 

4.2.2 Trajectory trapped around double-frequency orbit 

Below, we will show that the trajectory from the right panels of 
Fig.l involves free oscillation in addition to the two oscillations 
driven by the bars. The power spectrum for this trajectory is shown 
in the middle row of Fig. 4. It shows four major peaks. The two 
highest peaks at 0J2 = 195.9 km s -1 kpc -1 and Wi =231.7 km s -1 
kpc -1 are separated by twice the difference between the pattern 
speeds of the bars, 2(^2 — fii), which indicates that they occur at 
the frequencies of oscillations driven by the bars. Again, from i ll It 
one can obtain the underlying orbital frequency of the particle: Q — 
139.8 km s -1 kpc -1 , and from averaging the radial coordinate one 
can get an approximation for the guiding radius ro = 1.579 kpc. 
In the linear approximation, the orbital frequency at this guiding 
radius is 137.5 km s - kpc -1 , again consistent with the nonlinear 
solution. About 70% of the total power resides in these two peaks 
at frequencies of the driving by the two bars. 

Between the two highest peaks in the power spectrum, there 
are two other peaks of lower amplitude. They are located at fre- 



quencies k — 211.9 km s -1 kpc -1 and k' = 215.5 km s -1 kpc - , 
which are not commensurate with wi and Uz, and they contain 
about 20% of the total power. Note that ui\ — At' = At — u)i, which 
means that k and At' are not independent, and therefore that these 
two peaks indicate only one additional frequency on the trajectory. 
We postulate that this frequency is the free oscillation frequency. In 
the linear approximation of Sect. 3, this frequency corresponds to 
/to in G7J' which for the guiding radius ro = 1.579 kpc is equal to 
202 km s . This value is somehow different from « and «' mea- 
sured in the power spectrum, but it remains in their vicinity. More- 
over, qualitatively one should expect that at this guiding radius At 
and At' lie in between uii and u>2, because in the linear approxima- 
tion £1 — At/2 takes values between Qi and Q2 there. 

All other peaks in the spectrum occur at frequencies which are 
functions of U\, u>2 and At, and power contained in those peaks de- 
creases rapidly with frequency. Thus it appears that the orbit from 
the right panels of Fig. 1 is trapped around a double-frequency orbit 
and oscillates around it with the characteristic frequency At. 

4.2.3 Chaotic trajectory 

Finally, the power spectrum of the trajectory from the left panels of 
Fig.l is shown in the bottom rows of Fig.4. There are more peaks, 
and also much power contained in between the peaks. In order to 
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validate the significance of the peaks in this spectrum, we calcu- 
lated again power spectra for the three trajectories from Figs. 1 and 
2, but now followed for 600 alignments of the bars, with dynami- 
cal variables recorded at 2 16 equal time intervals. While the power 
spectra for the two trajectories analyzed above do not show signifi- 
cant changes, the new spectrum for the trajectory from the left pan- 
els of Fig. 1 , marked with the thick line in the bottom-left panel of 
Fig.4, shows multitude of peaks, each of them containing no more 
than 2% of the total power. Kandrup, Eckstein & Bradley (1997) 
showed that a multitude of frequencies indicates chaos (see also 
Athanassoula 2005). In the case analyzed here, there is a contin- 
uous distribution of power with broad maximum in the frequency 
range between 1 20 and 1 90 km s ~~ 1 kpc _ 1 . Such a power spectrum 
is characteristic of a chaotic trajectory. 

4.3 The vicinity of a double-frequency orbit 

In order to study how particles get trapped around double- 
frequency orbits, one can explore a range of initial conditions 
around the ones that characterize it. We developed an automatic 
procedure that finds local maxima (peaks) in the power spectrum at 
frequencies uj\ and U2 corresponding to the driving by the two bars 
and at the free oscillation frequency. We based this procedure on 
the fact that the two frequencies of driving by the two bars are al- 
ways separated by 2(^2 — fii ), where Q,\ and Q2 are the (constant) 
pattern speeds of the two bars. Thus, the algorithm first searches for 
the highest peak in the power spectrum, and records its frequency 
lo. Then it examines power spectrum at frequencies u±2(fi2 — fit)- 
If there is a peak in one of these frequencies, then this frequency 
and uj are recorded as lui and L02 on this orbit, as defined in (TTJ , 
Furthermore, the free oscillation frequency k is found by searching 
for the peak with highest power that does not occur at a frequency 
which is a linear combination of u>i and UJ2- 

In our analysis, we kept focus on a particle that starts at the 
1.7-kpc distance from the galaxy centre, on the minor axis of the 
aligned bars, with its initial velocity perpendicular to that axis. We 
analyzed a range of values of the initial velocity between 130 km 
s _1 and 230 km s _1 . The initial conditions for the three trajecto- 
ries studied above fall within this range. The dotted line in the upper 
panel of Fig. 5 indicates the average radius ro = (£r l )/JV (where 
the meaning of the symbols is the same as in Sect.4.2) as a func- 
tion of the initial velocity v x . As expected, higher initial velocity 
results in larger average radius, although the relation becomes not 
monotonic and noisy at initial velocities above 205 km s" 1 . When 
the average radius is taken as the radius of the guiding centre in 
the linear approximation, one can calculate the orbital frequency Q. 
and the free oscillation frequency k in this approximation. They are 
drawn in thick solid and dashed lines, respectively. 

On the other hand, once we derive wi and LU2 from the power 
spectrum of the trajectory, the orbital frequency Q. in the full non- 
linear solution can be obtained from i ll It . The values of Q derived 
from the measured values of Wi and L02 are the same, and they 
are plotted with a thin solid line in the upper panel of Fig.5. They 
closely follow the values derived in the linear approximation for 
most of the range, although there are clear departures at initial ve- 
locities above 205 km s" 1 . Although our algorithm to analyse the 
power spectrum allows us to extract automatically the free oscilla- 
tion frequency «, one should also expect a peak at k' — u>\ +0J2 — K, 
because two frequencies, wi and L02, are already present in the sys- 
tem. In fact, the algorithm finds peaks at two frequencies k and 
k', whose sum is wi + u>2. Both of these frequencies are drawn in 
the upper panel of Fig.5 (thin dashed lines), since from the power 
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Figure 5. Parameters that characterize orbits in Model 1 by Maciejewski & 
Sparke (2000), which start at 1.7-kpc distance from the galaxy's centre, on 
the minor axis of the bars when the bars are aligned, and with a velocity 
vector at right angle to it, plotted as a function of the magnitude of the stall- 
ing velocity v x . Upper panel: Thick dotted line marks the average radius 
of the orbit, ro = (Er t )/N. The epicyclic approximation for the angular 
frequency Q and for the frequency of free oscillations k at this radius are 
plotted with thick solid and dashed lines, respectively. The values of Q and 
k derived from the power spectrum of the orbit are plotted with thin solid 
and dashed lines, respectively. Lower panel: Percentage of the total power 
residing at the frequencies of the major peaks. The thin dashed line marks 
the power P(Q\) at the driving frequency of the outer bar, and the thin 
dotted line marks P(Q2) at the driving frequency of the inner bar. The sum 
of P(Qi) and P(Q2), ^(^1+2), is plotted with a thick dashed line. The 
thick dotted line marks the power P(k) in the two peaks arising from the 
free oscillation frequency. The thick solid line marks the median width of 
the ring enclosing particle positions at consecutive alignments of the bars, 
in units of the average radius of the orbit ro (see Sect.4.4). 



spectrum one cannot determine which one is the actual k. Their 
functional dependence on the initial velocity v x is similar to that 
obtained in the linear approximation, although the slope is slightly 
different, and the quantitative discrepancy is larger than in the case 
of the orbital frequency Q. Again, k and k' as a function of the 
initial velocity become ragged for v x > 205 km s _1 . 

One can also examine how much of the total power resides 
in various peaks in the Fourier spectrum. In the lower panel of 
Fig.5, a thin dashed line indicates the power P(Qi) residing in the 
peak at the frequency of driving by the outer bar, while a thin dot- 
ted line indicates the power P(fi2) in the peak at the frequency 
of driving by the inner bar. With the thick dashed line we plotted 
the power P(fli+2) residing in both peaks above. Power P(k), 
residing in the two peaks related to the free oscillation frequency, 
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is plotted with a thick dotted line. From the plot one can see that 
in the range of initial velocity v x between 130 km s _1 and 205 
km s _1 , P(fii+2) and P(k) sum up roughly to 100%. This means 
that only three frequencies are present on those trajectories, as ex- 
pected for the regular orbits (see Sect.2). In this range, the larger 
the value of P(k), the smaller P(fii+2). However, throughout this 
range P(f2i+2) > 0.5, which means that more than half of the to- 
tal power is located in the two peaks at the driving frequencies of 
the two bars. Orbits with larger average radii ro have more power 
in the driving frequency of the outer bar, while orbits with smaller 
ro have most of their power in the driving frequency of the inner 
bar. For v x around 180 km s _1 , up to 95% of the total power re- 
sides at the driving frequencies of the two bars, and only 0.03% 
in the free oscillations. Thus the trajectory with v x around 180 km 
s _1 is made up almost entirely of oscillations with the driving fre- 
quencies of the two bars. This is the double-frequency orbit from 
Fig.2. The more v x differs from 180 km s _1 , the larger fraction of 
the total power resides in the free oscillations: up to 38% and 33% 
at v x = 130 km s _1 and v x = 202 km s _1 , respectively. This is 
typical for trajectories trapped around stable regular orbits. An ex- 
ample of such a trajectory is the one from the right panels of Fig. 1 , 
for which v x = 150 km s _1 . 

For v x > 205 km s _1 , the summed power in peaks related 
to the three frequencies considered here is far smaller than 100%. 
The initial conditions of the trajectory from the left panels of Fig.l 
belong to this range of v x . The power spectrum of that trajectory 
(Fig.4, bottom row) shows the power distributed in a continuum of 
frequencies, which is characteristic for a chaotic orbit. Thus trajec- 
tories with v x > 205 km s _1 are most likely chaotic. This explains 
the noise for v x > 205 km s -1 in the curves that mark the frequen- 
cies derived from Fourier analysis in the top panel of Fig.5. One 
can point out the region of onset of chaos at v x between 200 km 
s _1 and 210 km s~\ where both P(f2r+2) and P(k) drop sharply 
with v x . 

4.4 Construction of loops from double-frequency orbits 

In the epicyclic approximation, the double-frequency orbits have a 
nice feature, namely that particles following them populate loops: 
closed curves that return to their original shape and position at ev- 
ery alignment of the bars. One may therefore expect that also in the 
general case considered here particles on double-frequency orbits 
populate loops. In order to check, whether it is the case, we fol- 
low the same procedure as in the epicyclic approximation. We start 
a particle from the minor axis of the aligned bars, and record its 
position and velocity every time the bars align again. Let xo mark 
the starting position and velocity of the particle, x\ its position and 
velocity at the first consecutive alignment of the bars and so on. 
Positions of a particle on a double-frequency orbit are overplotted 
in Fig.2 for the first 20 alignments of the bars. In the left panel of 
Fig.6, we plotted 60 positions for 60 alignments. They seem to be 
arranged on a closed curve of oval shape. 

Let us now take N values of initial conditions for a set of 
N particles at one alignment of the bars. By construction, a particle 
starting at Xi will be located at Xi+i at the next alignment of the 
bars. Thus, if positions of a particle at every alignment of the bars 
populate a closed curve, this curve will transform into itself every 
time the bars align. This is what we defined as a loop. In the left 
panel of Fig.6 one can see that positions of a particle on a double- 
frequency orbit, written at every alignment of the bars, populate a 
closed curve, which is a loop. Therefore, one may expect that the 
relation between double-frequency orbits and loops found in the 




Figure 6. Left panel: Positions of a particle on the double-frequency or- 
bit from Fig.2 at consecutive alignments of the bars generate a loop. Here 
60 positions are recorded. Remaining panels: Evolution of the loop during 
one relative period of the bars. The bars, outlined with solid lines, rotate 
counterclockwise. 



linear approximation also holds in the general case. Note that this 
construction of a loop requires the particle to be followed for many 
relative periods of the bars, hence it will likely produce loops that 
are maps of stable double-frequency orbits only. 

The loop returns to its original shape and position at every 
alignment of the bars. Recording positions of particles from the 
loop at any time between alignments of the bars will unveil how 
the loop evolves between the alignments. This has been shown in 
the remaining panels of Fig.6, where we plotted locations of these 
particles at three other relative positions of the bars. The particles 
remain confined to the loop, but the shape of the loop varies in 
between the alignments. 

Will particles on other orbits, which involve free oscillations, 
take any particular positions at every alignment of the bars? One 
may expect that since such orbits oscillate around double-frequency 
orbits, the recorded positions will gather around the loop. If so, they 
should be enclosed within a band of certain width that increases 
with the amplitude of free oscillation. This band, surrounding the 
loop, is closed and, following the notation introduced by Maciejew- 
ski & Sparke (2000), we will call it the ring. In order to see how 
the width of the ring changes with the amplitude of free oscillation, 
we applied the following algorithm, similar to the one used by Ma- 
ciejewski & Sparke (2000). We registered particle positions on the 
orbit at 400 alignments of the bars. We divided the full 2-7T angle 
into 80 equal parts, and for each part we calculated the range Ar of 
radii among the positions falling into this part. The median of the 
Ar values is taken as the estimate of the ring width W. 

In the bottom panel of Fig.5 we plot the ring width W, nor- 
malized to the average radius of the orbit ro, as a function of the 
starting velocity. At velocities between 130 km s _1 and 200 km 
s _1 , the ring width W is below 30% of the average radius, hence 
recorded positions occupy a well defined ring. The ring gets thin- 
ner for velocities close to 180 km s _1 . This is the initial condition 
for the double-frequency orbit from Fig.2 that generates the loop 
from Fig.6. For this orbit W is about 2%. Generally, in the range 
of initial velocities considered here, W behaves in a very similar 
way to P(k). The larger the free oscillations around the double- 
frequency orbit, the wider the ring enclosing recorded positions. 
Therefore, we confirm our expectations that orbits trapped around 
double-frequency orbits map onto rings surrounding the loop. 

Note that in this algorithm, W is zero only for circular loops or 
under infinitely fine sampling in angle. Otherwise, W is only near 
zero for a loop. However, unlike the amplitude of free oscillation 
with frequency k, W is very straightforward to measure and to vi- 
sualize. The excellent agreement between the functional behaviour 



10 Witold Maciejewski and E. Athanassoula 



of W and P(n) exemplified in Fig.5 makes W a useful indicator 
of a loop. 

For initial velocities between 200 km s _1 and 210 km s _1 the 
median ring width sharply rises from 30% to 70% of the average 
radius. Such large median width indicates that recorded positions 
are no longer confined by a ring, but rather densely populate some 
fraction of the plane, and do not gather on any curve. This is char- 
acteristic for chaotic orbits, which is consistent with the frequency 
analysis above. 



5 DISCUSSION AND CONCLUSIONS 

The goal of this paper is to present an argument that the backbone 
of an oscillating potential of a doubly barred galaxy is built out of 
double-frequency orbits. In Sect.3, by using the linear approxima- 
tion, we showed that two frequencies are indispensable in motion 
in such systems. A single frequency, and thus a closed periodic 
orbit, can occur only when these two frequencies are commensu- 
rate. However, even in such case they are two separate frequencies, 
which happen to have the same numerical value. This is a condition 
for resonance, like the one in the single bar when the frequency of 
free oscillations is equal to that of driving by the bar. 

Double-frequency orbits in the oscillating potential of double 
bars do not close. Their appearance is also different in the reference 
frame of each bar. However, their perception and analysis is facili- 
tated greatly by the fact that they map onto loops, which are closed 
curves. The mapping is done by recording the positions on the orbit 
at every alignment of the bars. In such a construction, points that 
are on the loop will remain on the loop. A loop can be viewed as 
a set of particles, and therefore its appearance is independent from 
the frame of reference. However, the loop oscillates in time with a 
period equal to the relative period of the bars. 

In the linear approximation we showed rigorously that each 
double-frequency orbit maps onto a loop. In the general, non-linear 
case, we found that when the Fourier analysis of the orbit indicates 
almost all power at the frequencies of driving by the bars, and al- 
most no power at the the free oscillation frequency, the map of the 
orbit is to a good approximation a 1-D curve. This is the loop. When 
free oscillations are present, the orbit maps onto a set of points in 
a ring surrounding the loop. The width of this ring is higher when 
more power is in the free oscillations. 

From the power spectrum of an orbit in the oscillating poten- 
tial of a doubly barred galaxy we can retrieve its three fundamental 
frequencies (frequencies of forcing by the two bars and the fre- 
quency of free oscillations) when the orbit is regular. Otherwise a 
continuous power spectrum is recovered when the orbit is chaotic. 
The Fourier algorithm that we apply displays the baseline at no 
more than 10~ 4 of the total power, hence peaks higher than that 
can be easily extracted. Thus, even if the power in e.g. the free 
oscillation is very small around the loop, our algorithm recovers 
its frequency correctly, as the continuity of the k(v x ) curve in the 
top panel of Fig.5 indicates. In our algorithm, the power spectrum 
of each trajectory is analyzed independently from power spectra 
of neighbouring trajectories, and frequencies recovered on one tra- 
jectory do not enter as an 'initial guess' for those on the nearby 
trajectory. Therefore, continuity of Q and k as functions of the ini- 
tial velocity v x is a genuine property of the system, and not of the 
algorithm. 

Methods employed in this paper in the search for double- 
frequency orbits and loops require the particle to be followed for 
many (i.e. hundreds) relative periods of the bars. We do not perform 



fine tuning of the initial conditions on the trajectory, hence recover- 
ing unstable orbits with these methods is rather unlikely. Therefore 
double-frequency orbits recovered with the methods employed in 
this paper are stable in the sense that they trap regular orbits, which 
wind tightly around them and map onto rings enclosing loops. 

We conclude that in a potential of two independently rotating 
bars, regular trajectories are trapped around parent stable double- 
frequency orbits. They oscillate around these orbits, like trajecto- 
ries in a single bar oscillate around stable closed periodic orbits. 
The two frequencies present on the parent orbits correspond to the 
frequencies of driving by the two bars. Therefore the parent orbits 
do not close in any reference frame. However, they can be mapped 
onto the loops, whose appearance is independent from the refer- 
ence frame. As such, loops are a unique tool that greatly facilitates 
the search for regular orbits in oscillating potentials. 
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